## [1] "ATTRIBUTE"  "WETLAND_TY" "ACRES"      "Shape_Leng" "Shape_Area"

(Stewart et al. 2024)

## # A tibble: 9 × 2
##   hpu   result
##   <fct>  <dbl>
## 1 Bh    0.290 
## 2 Bhs   0.118 
## 3 Bhsm  0.168 
## 4 Bi    0.227 
## 5 E     0.149 
## 6 H     0.583 
## 7 I     0.259 
## 8 O     0.0601
## 9 T     0.127
## Warning: Groups with fewer than two data points have been dropped.
## Warning: Removed 11 rows containing non-finite values (`stat_ydensity()`).
## Warning: Groups with fewer than two data points have been dropped.
## Warning: Removed 11 rows containing missing values (`geom_point()`).

Is there a model?

## 
## Call:
## lm(formula = WET ~ hpu, data = hb_pedons_ext)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -0.28180 -0.12423 -0.06280  0.06577  0.71024 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  0.28980    0.02539  11.413  < 2e-16 ***
## hpuBhs      -0.17225    0.04077  -4.225 3.08e-05 ***
## hpuBhsm     -0.12180    0.07571  -1.609  0.10858    
## hpuBi       -0.06304    0.03787  -1.665  0.09689 .  
## hpuE        -0.14123    0.04905  -2.879  0.00424 ** 
## hpuH         0.29320    0.21547   1.361  0.17450    
## hpuI        -0.03110    0.05777  -0.538  0.59076    
## hpuO        -0.22975    0.05647  -4.069 5.88e-05 ***
## hpuT        -0.16275    0.03300  -4.931 1.28e-06 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.214 on 339 degrees of freedom
##   (11 observations deleted due to missingness)
## Multiple R-squared:  0.1141, Adjusted R-squared:  0.09317 
## F-statistic: 5.456 on 8 and 339 DF,  p-value: 1.741e-06
## Analysis of Variance Table
## 
## Response: WET
##            Df  Sum Sq  Mean Sq F value    Pr(>F)    
## hpu         8  1.9984 0.249801  5.4562 1.741e-06 ***
## Residuals 339 15.5203 0.045783                      
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Warning: Removed 5 rows containing non-finite values (`stat_ydensity()`).
## Warning: Removed 5 rows containing missing values (`geom_point()`).
## `geom_smooth()` using formula = 'y ~ x'
## Warning: Removed 30 rows containing non-finite values (`stat_smooth()`).
## Warning: Removed 30 rows containing missing values (`geom_point()`).

Is there a model in W3?

## 
## Call:
## lm(formula = WET ~ hpu, data = w3_pedons_ext)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -0.15737 -0.08136  0.00100  0.05010  0.22100 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         0.23275    0.03892   5.981 2.59e-06 ***
## hpuBhs podzol      -0.06037    0.05504  -1.097   0.2827    
## hpuBimodal podzol  -0.14325    0.08702  -1.646   0.1118    
## hpuE podzol        -0.08355    0.06275  -1.331   0.1946    
## hpuLithic Histosol  0.03725    0.06740   0.553   0.5852    
## hpuTypical podzol  -0.11575    0.06275  -1.845   0.0765 .  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.1101 on 26 degrees of freedom
##   (5 observations deleted due to missingness)
## Multiple R-squared:  0.2335, Adjusted R-squared:  0.08605 
## F-statistic: 1.584 on 5 and 26 DF,  p-value: 0.1997
## Analysis of Variance Table
## 
## Response: WET
##           Df   Sum Sq  Mean Sq F value Pr(>F)
## hpu        5 0.095937 0.019187  1.5837 0.1997
## Residuals 26 0.315003 0.012115
## 
## Call:
## lm(formula = WET ~ water_table, data = w3_pedons_ext)
## 
## Residuals:
##        3        6       11       17       19       28       35 
## -0.20724 -0.11312  0.04029  0.06486  0.19226 -0.04112  0.06408 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)  
## (Intercept)  0.350684   0.158615   2.211    0.078 .
## water_table -0.002206   0.002690  -0.820    0.450  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.1445 on 5 degrees of freedom
##   (30 observations deleted due to missingness)
## Multiple R-squared:  0.1185, Adjusted R-squared:  -0.05776 
## F-statistic: 0.6724 on 1 and 5 DF,  p-value: 0.4495
## Analysis of Variance Table
## 
## Response: WET
##             Df   Sum Sq  Mean Sq F value Pr(>F)
## water_table  1 0.014031 0.014030  0.6724 0.4495
## Residuals    5 0.104337 0.020867

Watershed 3 wells

## Warning: package 'stringr' was built under R version 4.3.1

Other soil moisture data

##    plot         WET
## 1   IL1 0.111505258
## 2   IL2 0.017028480
## 3   IL3 0.016389744
## 4   IH1 0.014945904
## 5   IH2 0.300747536
## 6   IH3 0.009245352
## 7   E01 0.108972447
## 8   E02 0.043977207
## 9   E03 0.042729688
## 10  E07 0.010488627
## 11  E08 0.022509384
## 12  E09 0.148845048
## 13  E10 0.342572154
## 14  E12 0.022831787

## `geom_smooth()` using method = 'gam' and formula = 'y ~ s(x, bs = "cs")'
## Warning: Removed 7674 rows containing non-finite values (`stat_smooth()`).
## `geom_smooth()` using method = 'gam' and formula = 'y ~ s(x, bs = "cs")'
## Warning: Removed 19943 rows containing non-finite values (`stat_smooth()`).
## `geom_smooth()` using method = 'gam' and formula = 'y ~ s(x, bs = "cs")'
## Warning: Removed 32704 rows containing non-finite values (`stat_smooth()`).
## `geom_smooth()` using method = 'gam' and formula = 'y ~ s(x, bs = "cs")'
## Warning: Removed 15780 rows containing non-finite values (`stat_smooth()`).
## Warning: Removed 32704 rows containing missing values (`geom_point()`).
## Warning: Removed 19943 rows containing missing values (`geom_point()`).
## Warning: Removed 7674 rows containing missing values (`geom_point()`).
## Warning: Removed 15780 rows containing missing values (`geom_point()`).

Old Denitrification soil moisture & oxygen

##   ID       UPL        WET
## 1  1 0.9856686 0.01433141
## 2  2 0.9671750 0.03282499
## 3  3 0.9192316 0.08076836

## Warning: Removed 1 rows containing missing values (`geom_point()`).
## Removed 1 rows containing missing values (`geom_point()`).
## Removed 1 rows containing missing values (`geom_point()`).
## Warning: Removed 26 rows containing missing values (`geom_point()`).
## Warning: Removed 1 rows containing missing values (`geom_point()`).
## Warning: Removed 90 rows containing missing values (`geom_point()`).
## Warning: Removed 1 rows containing missing values (`geom_point()`).
## Removed 1 rows containing missing values (`geom_point()`).
## Removed 1 rows containing missing values (`geom_point()`).

## Warning: Removed 209 rows containing missing values (`geom_point()`).
## Warning: Removed 74 rows containing missing values (`geom_point()`).

## Warning: Removed 4989 rows containing missing values (`geom_point()`).

## `geom_smooth()` using formula = 'y ~ x'
## `geom_smooth()` using formula = 'y ~ x'

Stewart, Anthony J., Meghan Halabisky, Chad Babcock, David E. Butman, David V. D’Amore, and L. Monika Moskal. 2024. “Revealing the Hidden Carbon in Forested Wetland Soils.” Nature Communications 15 (11): 726. https://doi.org/10.1038/s41467-024-44888-x.